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ABSTRACT 


The explosive shock created by the underwater explosion of a mine or torpedo in 
close proximity to a surface ship can severely threaten the combat capability and 
survivability of the ship. MIL-S-901D specifies the shock test procedures and acceptance 
criteria for all shipboard systems that must resist high impact mechanical shock. While the 
U.S. Navy’s Mediumweight Shock Machine with its standard equipment mounting fixture 
can subject a combat systems component to more severe shock excitations than 
experienced in actual ship shock trials, it cannot simulate the lower frequency excitations 
typically transmitted through a ship’s superstructure during shock trials that expose 
equipment to catastrophic resonant vibration. This study is an experimental investigation 
into the dynamic response of the recently built Two Degree-of-Freedom (2DOF) Tuned 
Deck Simulator (TDS) for the Mediumweight Shock Machine (MWSM) to evaluate its 
potential role in the pre-acceptance shock qualification of new shipboard combat systems 
equipment. Upon completion of final characterization testing, the 2DOF-TDS could be 
integrated into the mediumweight shock qualification procedures of MIL-S-901D. This 
improvement could significantly enhance the capacity of a warship to absorb damage and 


still maintain its mission integrity. 
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I. INTRODUCTION 


Survivability is defined as the capacity of a ship to absorb damage and maintain 
mission integrity. In September 1988, the Chief of Naval Operations issued a 
comprehensive policy directive to incorporate survivability features into new surface ship 
designs, ship overhauls, and new combat systems. This directive ranked survivability 
equal in significance to such fundamental design requirements as maneuverability and 
combat systems capability. The minimum baseline of ship survivability includes: armor, 
shielding and signature reduction, and the shock hardening of vital systems. The directive 
tasked the Commander, Naval Sea Systems Command, with establishing and validating 
survivability performance and test standards matched to actual and anticipated threats. 
(OPNAV, 1988, pp. 2-5) 

One major threat to surface ships is the explosive shock created by the underwater 
explosion of a mine or torpedo in close proximity to the ship. The U.S. Navy’s shock 
hardening policies and procedures are primarily concerned with ensuring the operability of 
mission essential systems. Shock hardening is required for all ships designated as capable 
of operating in the combat shock environment, including: surface combatants, aircraft 
carriers, amphibious ships, mine warfare ships, and others designated to operate inside 
combat zones. The Naval Sea Systems Command (NAVSEA) requires that all new 
construction ships be shock hardened and that all overhaul contracts and major combat 
system design upgrades include appropriate shock qualification requirements. (NAVSEA, 


1989, p. 4) 


The U.S. Navy’s shock qualification requirements are contained in Military 
Specifications, MIL-S-901D, “Shock Tests, High Impact; Shipboard Machinery, 
Equipment and Systems, Requirements for,” which specifies the shock test procedures and 
acceptance criteria for all grades and classes of shipboard systems that must resist high 
impact mechanical shock. The three methods of shock testing are: 

e Lightweight Shock Machine Test for equipment and mounting with a 

combined weight less than 550 pounds using the Lightweight Shock Machine 
(LWSM) 

e Mediumweight Shock Machine Test for equipment and mounting with a 

combined weight less than 7400 pounds using the MWSM 

° Heavyweight Shock Test for equipment weights exceeding MWSM capacity 

The Lightweight Test requires the LWSM which consists of a welded steel 
framework equipped with a rotating anvil plate and two 400 pound hammers that can be 
configured for shock testing equipment along three orthogonal axes without remounting 
(Clements, 1972, p. 32). The Heavyweight Test requires a floating shock platform and 
submerged explosive charges and is recommended in cases where the simulation of deck- 
mounted conditions is desired. (MIL-S-901D, 1989, pp. 12-16) 

This research is focused on the Mediumweight Shock Machine (MWSM) and the 
new Two Degree-of-Freedom Tuned Deck Simulator (2DOF-TDS) upon which a test 
item would be base-mounted for shock loadings in the vertical direction. Before 
proposing a future role for the 2DOF-TDS in the shock qualification of vital combat 


systems, a description of current mediumweight shock qualification is necessary. 


A schematic drawing of the MWSM is provided in Figure 1. The MWSM is a 
hammer-anvil table apparatus originally designed in 1942. After release from a specified 
height, the 3000 pound hammer swings through an arc of up to 270 degrees and strikes 
the 4500 pound anvil table from below, imparting an upward acceleration to both the table 
and the test item mounted on the table. The anvil table vertical travel distance after impact 
is limited to either 3.0 inches (due to 12 foundation bolts which act as mechanical stops), 
or 1.5 inches (if installed pneumatic jacks are used to raise the anvil table prior to impact). 


(Clements, 1972, p. 69) 


All items subject to mediumweight shock testing receive a minimum of six hammer 
blows which are applied in three groups of two blows each. MIL-S-901D specifies the 
height of the hammer drop and the anvil table travel distance for each group. One blow 
of each group is conducted with the test item in its normal attitude to simulate shipboard 
shock loadings in the vertical direction. The second blow of each group is conducted with 
the test item mounted in the inclined position to simulate shipboard shock loadings in the 
athwartship direction. Test items that are essential to the safety and combat capability of 
the ship must be operated in their normal operating modes atop the MWSM during the 
first and third groups of blows. Alternate modes of operation (i.e., “standby” mode), are 
tested during the second group of blows. (MIL-S-901D, 1989, p. 19) 

According to MIL-S-901D, base-mounted items intended for horizontal surfaces 
aboard ship are normally mounted on the anvil table by means of a standard fixture, 
consisting of a specified number of support channels atop two shipbuilding rails attached 


to the anvil table, for shock loadings in the vertical direction. Figure 13 of MIL-S-901D 


illustrates typical standard fixture arrangements and provides a table which enables the 
user to determine the number of supporting fo::: inch car building channels required for a 
given equipment weight and size (MIL-S-901D, 1989, pp.67-69). A standard inclined 
fixture for shock loading in the athwartship direction is used to incline the test item in the 
same direction that the item would incline aboard ship if the ship were to roll 30 degrees. 

All standard mounting fixtures are intended to represent hull-mounted conditions 
that exist at main structural members of the ship (primarily below the main deck). For 
deck-mounted conditions, MIL-S-901D with some exceptions limits the use of standard 
mounting fixtures to equipment installed aboard ship without resilient mountings (ie., 
Class I equipment). primarily located on the main deck and above, but also on decks and 
non-structural bulkheads below the main deck. When simulating deck-mounted 
conditions, MIL-S-901D stipulates that the fundamental frequency (with the equipment 
installed on the fixture) shall not be lower than 25 Hz in each principal shock direction. 
(MIL-S-901D, 1989, pp. 2-12) 

The primary issue is whether or not the MWSM and standard fixture can realistically 
simulate a ship’s actual response to an underwater explosion for the pre-acceptance shock 
testing of new shipboard combat systems, particularly equipment designated for deck- 
mounting in a ship’s superstructure. While the MWSM and standard fixture generate a 
more severe shock excitation than experienced in actual ship shock trials, the test item is 
not subjected to the lower frequency excitations typically transmitted through a ship’s 
superstructure during shock trials (Cox, 1993, pp. 7-8). Ifthe ship’s structure excites the 


combat systems equipment at lower frequencies than those generated by the MWSM and 


standard fixture, an actual underwater explosion in close proximity to the ship could 
expose the equipment to catastrophic resonant vibration not previously simulated during 
pre-installation shock qualification tests. A 1992 study revealed that the short-duration, 
high-impulse, high-frequency excitation waveform generated by the MWSM_ and 
transmitted to the test item via the standard fixture significantly differs from actual 
waveforms observed at various equipment locations aboard ship during shock trials 
(Corbell, 1992, p. 4). For more accurate mediumweight shock testing of vital combat 
systems equipment, Corbell proposed that the standard fixture be replaced by a two 
degree-of-freedom mounting fixture pre-tuned to simulate the dominant shock 
characteristics derived in the ship class pre-shock trial analyses. 

The 2DOF-TDS for the MWSM was designed by Cox, and a prototype was built in 
1993. This study is an experimental investigation into the dynamic response of the 2DOF 
Tuned Deck Simulator to evaluate its potential role in mediumweight shock qualification. 
Chapter II provides background information on the design and construction of the 2DOF- 
TDS. The characterization of the 2DOF-TDS on the MWSM through actual shock 
testing is presented in Chapter III. Chapter IV addresses the possibility of incorporating 
the 2DOF-TDS into mediumweight shock qualification procedures as an addendum to 
MIL-S-901D, The conclusion and recommendations toward addendum validation are 


presented in Chapter V. 


Il. BACKGROUND 


Foremost among the U.S. Navy’s surface ship survivability priorities is the shock 
hardening of mission essential equipment. Given the increasing budgetary constraints 
facing the U.S. Navy today, making incremental improvements to the simulation accuracy 
of the six government-owned MWSM’s may be more feasible than developing a 
replacement for the MWSM. Before exploring the design and construction of the 2DOF 
Tuned Deck Simulator, a comparison of the DDG-51 Class pre-shock trial spectra with 


that generated by the MWSM and standard fixture is in order. 


A. DDG-51 PRE-SHOCK TRIAL SPECTRA VERSUS 
MWSM SPECTRA 


In underwater shock analysis, the shock spectrum is defined as the maximum 
absolute response of an undamped single DOF system subjected to a base acceleration 
initiated by contact with a shock wave. The shock spectra encompass the undamped 
single DOF system’s absolute maximum response envelope over a wide range of system 
natural frequencies resulting from a specific excitation. For a given excitation, the 
resultant shock spectra reveal the peak resonance responses experienced by the system. If 
the complex combination of natural frequencies in a ship’s structure following an 
underwater explosion further excites the natural frequencies of components inside vital 
combat systems equipment, significant equipment damage and a drastic reduction in the 


ship’s combat capability may result. In order to simulate this shipboard excitation on the 


MWSM, it is not necessary to reproduce the actual waveform, only its frequency 
components at their respective acceleration amplitudes (Cox, 1993, p. 9). 


1.‘ Transient Shock Analysis of DDG-51 Class Forward Deckhouse 


In 1992, Costanzo and Murray of the David Taylor Research Center 
Underwater Explosion Research Division (DTRC/UERD), completed a transient shock 
analysis of the DDG-51 Class forward deckhouse and mast in preparation for the shock 
trial of DDG-53 in 1994. Their finite element model of the forward deckhouse included 
three lumped, rigid masses representing significant combat systems components, ranging 
in weight from 325 to 4600 pounds, all located on the ship’s 0-3 level. Applying the 
maximum shock trial excitations obtained during previous CG-47 Class cruiser shock trials 
to major support bulkheads of the DDG-51 forward deckhouse model, Costanzo and 
Murray computed the dynamic responses, including the predicted vertical acceleration 
waveforms and the resultant shock spectra, at the approximate locations of the three 
combat systems components. (Corbell, 1992, pp. 11-15) 

Although ail three components were located on the 0-3 level of the DDG-51 
Class forward deckhouse, the predicted acceleration waveforms and their respective shock 
spectra were dramatically different. Corbell compared the Costanzo and Murray finite 
element model’s acceleration waveform with that generated by Clement’s MWSM model 
and concluded that the MWSM with standard fixture could not reproduce the 
“characteristic shock phenomena observed in the field” (Corbell, 1992, pp. 15-16). 

2. Need for an Improved MWSM Fixture 
To substantiate Corbell’s conclusion, a comparison between Costanzo and 


Murray’s predicted acceleration time history for one of their three tested components, the 
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SPY-1D Radar Beam Programmer weighing 1000 pounds, and actual MWSM (with 
standard fixture) acceleration data obtained by the Naval Undersea Warfare Center 
(NUWC), New London, in a June 1993 test series to characterize the MWSM was 
performed (depicted in Figure 2). The NUWC test data chosen for this comparison 
involved a 1000 pound test weight supported by a 300 pound aluminum plate atop the 
standard fixture (Messina, 1993). Although the combined weight on the MWSM 
exceeded the weight of the beam programmer by 300 pounds, the difference in the 
respective acceleration time histories is significant. 

More significant than this waveform comparison in the time domain is a 
comparison of the shock wave characteristics displayed in the frequency domain as shock 
spectra. The respective absolute acceleration shock spectra are depicted in Figure 3. 
Although some correlation of the respective shock spectral peak locations exists, the 
acceleration amplitudes in the MWSM’s shock spectra are significantly higher than those 
of the beam programmer. While MIL-S-901D does not require that any particular 
waveform or spectra be reproduced during shock testing, this comparison clearly signifies 
the need for an improved MWSM fixture that can more accurately simulate shock 
phenomena observed in the field for the pre-acceptance shock qualification of new 
combat systems equipment. 

From the DTRC/UERD DDG-51 Class pre-shock trial analysis, Corbell 
observed that the vertical orientation shock spectra of both light and medium weight 
equipment revealed two common, well-defined dominant peaks at 60 Hz and 155 Hz, 


respectively. Corbell mathematically modeled a 2DOF uniaxial mounting fixture that was 


tuned to achieve 60 Hz and 155 Hz as its coupled response frequencies and concluded that 


the design and construction of such a fixture was warranted. (Corbell, 1992, pp.51-73) 


B. DESIGN AND CONSTRUCTION OF THE 2DOF 
TUNED DECK SIMULATOR 


1. Initial Design Phase 

Following Corbell’s lead, Cox set out to design a tunable 2DOF Deck 
Simulator for the MWSM that would excite combat systems equipment with most of the 
energy concentrated at two dominant resonant frequencies pre-selected from actual ship 
shock trial spectral analysis or from modal analysis at the equipment’s designated location 
aboard ship. Initially using Corbell’s target frequencies of 60 Hz and 155 Hz for shock 
spectral peaks, Cox modeled the 2DOF Tuned Deck Simulator as two mass-spring- 
damper systems coupled together as depicted in Figure 4. The mass “M1” represented the 
combined masses of the test object and associated mounting hardware. The mass “M2” 
represented the combined masses of the intermediate structure and associated hardware. 
The characteristic stiffness and damping properties of the two tiers were designated by 
“K” and “C” respectively. Cox later included the effective spring mass of his design’s 
upper spring beams into M1 and that of both upper and lower spring beams into M2. Cox 
also modeled the anvil acceleration associated with a hammer release height of one foot as 
a half-sine wave vertical acceleration impulse of 1 millisecond duration as done previously 


by Clements and Corbell. (Cox, 1993, p. 16) 


2. Model Refinement Phase 

After deriving the model’s coupled differential equations of motion in matrix 
form with coordinates relative to the motion of the anvil table, Cox calculated the first two 
natural frequencies. Defining f; and /> as the “uncoupled” tier natural frequencies, and /;, 
and f,2 as the “coupled” system natural frequencies, Cox derived equations relating the 
system’s coupled and uncoupled natural frequencies. Choosing 60 Hz and 155 Hz as his 
desired system natural frequencies, Cox iteratively solved these equations for the required 
uncoupled natural frequencies. For example, assuming a mass ratio (i.e., Mla), of 1.0 
and an uncoupled lower tier natural frequency (f2) of 100 Hz, an upper tier natural 
frequency (/;) of 94 Hz would be required to provide system natural frequencies (f,, and 
Jn2) of 60 Hz and 155 Hz. (Cox, 1993, pp. 18-20) 

Cox then uncoupled the model’s equations of motion, assumed values of 
damping between 0.03 and 0.05, and - using a MATLAB?’ simulation algorithm — obtained 
system displacement, velocity and acceleration responses to the half-sine wave impulse. 
Cox used this algorithm to determine the displacements of each tier relative to the anvil 
table for various MWSM hammer heights and weights atop the 2DOF-TDS. From this 
data, the maximum dynamic relative displacement within the 2DOF-TDS was obtained for 
a Static stress analysis of the components to aid in their material selection. (Cox, 1993, 
pp. 23-28) 


3. Final Design Phase and Stress Analysis 


With the objective of defining a structure that would both provide the desired 


response without plastic yielding or cracking during shock tests and include some of the 
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MWSM’s current accessories to minimize fabrication cost, Cox arrived at the two tier 
beam spring design pictured in Figure 5. Simply-supported I-beams with two equal 
concentrated forces symmetrically placed were chosen for the 2DOF Tuned Deck 
Simulator’s spring elements to achieve the required stiffness while minimizing weight. 
The existing shipbuilding rails atop the anvil table were chosen to support the four lower 
tier spring beams. An intermediate structure was placed between the three upper tier 
spring beams and the four lower tier spring beams to achieve the required uncoupled tier 
frequencies and mass ratio that define the system response, and to prevent unwanted 
deflections and vibrations from interfering with the desired response. The intermediate 
structure was designed with the same mounting geometry as the shipbuilding rails on the 
anvil table to enable the upper and lower spring beams to be interchanged. The test object 
was mounted on two “equipment specific” I-beams clamped across the upper spring 
beams. (Cox, 1993, pp. 32-38) 

Cox focused the stress analysis of the 2DOF Tuned Deck Simulator on the 
bending stress in the spring beams and clamps and on the tensile stress in the fasteners. 
Using his simulation algorithm of the 2DOF-TDS, loaded to capacity and excited with 
Clement’s maximum anvil table acceleration of 200 g’s (corresponding to a five foot 
hammer drop height), Cox obtained the maximum relative displacement of the spring 
beams and calculated the maximum bending stress in the spring beams to be 47 ksi (13 ksi 
in the clamps). The maximum tensile stress in the fixture bolts was 7.5 ksi. (Cox, 1993, 


p. 43) 
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Cox’s design drawings of the 2DOF Tuned Deck Simulator were approved by 
NAVSEA and forwarded to the Naval Surface Warfare Center, Silver Springs, Md. for 
construction. The bill of material for the 2DOF-TDS and a schematic drawing are 
included in Figure 6. 

4. Tunability Features of the 2DOF Tuned Deck Simulator 

Adding tunability features to the 2DOF-TDS design entailed providing the 
means to obtain specific system natural frequencies for a variety of equipment weights by 
altering the stiffness-to-mass ratio of the upper and lower tiers. The upper tier frequency 
can be adjusted by removing the center spring beam from the upper tier and/or by altering 
the top load spacing (of the equipment support I-beams) from each end of the upper 
spring beams. The lower tier frequency is adjusted by adding ballast to the intermediate 
structure as necessary (without exceeding the anvil table’s load limit of 7400 pounds). 
Cox utilized a MATHCAD® worksheet to determine the achievable system natural 
frequencies for various upper tier weights and 2DOF-TDS configurations. (Cox, 1993, 


pp. 36-47) 


C. DDG-51 MODAL TESTING AND TARGET 
RESONANT FREQUENCIES FOR THE 2DOF 
TUNED DECK SIMULATOR 


In September 1993, the Bath Iron Works Corporation conducted modal testing of 
the decks inside the combat systems equipment rooms of three DDG-51 Class ships to 


identify the lowest natural frequency of each deck. The transmissibility modal tests, in 


which accelerations at respective deck grid points were referenced to accelerations at a 
fixed point on each deck, revealed deck fundamental natural frequencies between 16.5 Hz 
and 25 Hz (Horton, 1993, p. 20). Better quality data was later obtained during a 
frequency response function modal test in Combat Systems Equipment Room 3 (CSER3), 
where accelerations at deck grid points were referenced to the input force applied to the 
deck. Using STAR® Modal Analysis algorithms, Horton identified four natural 
frequencies between 14.8 Hz and 29.6 Hz and their mode shapes in the deck of CSER3 
(Horton, 1993, p. 27). 

Horton’s findings and the existence of the coil spring Soft Deck Simulator, used for 
mediumweight shock testing equipment deck-mounted in submarines, with a natural 
frequency range of 20 to 25 Hz, led Cox to alter his target resonant frequencies for the 
2DOF Tuned Deck Simulator from 60 Hz and 155 Hz to 30 Hz and 80 Hz. Cox also 
introduced a semi-definite 3DOF model that accounted for anvil table interaction with the 
deck simulator approximately 50 milliseconds after hammer impact due to the sharp “table 
reversal” experienced by the anvil table when meeting its vertical travel limit. The 3DOF 
model’s shock spectra for various tier weights and a simulated hammer height of five feet 
displayed higher system natural frequencies at lower acceleration amplitudes than that of 
the 2DOF model with the same initial conditions. Cox concluded that the actual system 
natural frequencies should fall somewhere between those of the 2DOF and 3DOF models, 
respectively (Cox, 1993, p. 63). 

The 2DOF-TDS was built and delivered to the Naval Undersea Warfare Center, 


Newport, R.I. in September 1993. The first phase of fixture characterization testing on 
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the MWSM began on 19 October 1993. A photograph of the 2DOF-TDS with the 
aluminum mounting plate and a 500 pound test weight atop the MWSM anvil table is 
provided in Figure 7. The test procedure, test results and a comparison between the 


2DOF-TDS and the standard fixture in terms of shock spectra follow in Chapter III. 
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Til. CHARACTERIZATION OF THE 2DOF TUNED 
DECK SIMULATOR 


A. DESCRIPTION OF 2DOF TUNED DECK 
SIMULATOR TESTING 


The experimental investigation into the use of the newly fabricated 2DOF-TDS 
involved two test phases, both conducted at NUWC, Newport, R.I. The objective of the 
first test phase (19-20 October 1993) was to collect data on the 2DOF-TDS during shock 
testing on the MWSM with test loads of 500, 1000 and 1500 pounds. These test loads 
atop the fixture were chosen to match load conditions of previous MWSM testing that 
characterized the MWSM standard fixture in June 1993 (Messina, 1993). In accordance 
with the MWSM Test Schedule (Table 1 in MIL-S-901D), each of the six tests in Test 
Phase I consisted of three groups with two blows, or hammer drops, per group. 

The objective of the second test phase (01-02 February 1994) was to increase the 
accuracy of the acceleration shock spectra for six of the original 18 groups of blows by 
conducting six tests of ten blows each and averaging the respective acceleration shock 
spectra data. The six tests of Phase II were selected based on how well-defined the 
respective acceleration shock spectra peaks appeared in Phase I test results and on how 
well the Phase II results would best represent the diversity of hammer heights and anvil 
table travel limits specified in Table 1 of MIL-S-901D for the three different test weights 


atop the 2DOF-TDS. 
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The characterization of the 2DOF Tuned Deck Simulator was based on the 
acceleration shock spectra obtained in Test Phase II. For this reason, the description of 
the test procedure, presentation of test results and the comparison between the 2DOF- 
TDS and the standard fixture which follow will consider only Phase II data. 


1. Test Schedule 


The first three tests of Phase I] were conducted with the 2DOF-TDS in its 
“3/4” configuration, meaning all three upper spring beams and four lower spring beams 
were in place. The latter three tests were conducted with the 2DOF-TDS in its “2/4” 
configuration, meaning the center spring beam in the upper tier was removed while the 
lower tier was unchanged. 

The actual test schedule for Phase II is provided in Table 1 of Appendix B. 
Ten blows per test were conducted at hammer heights and anvil travel limits specified in 
Table 1 of MIL-S-901D based on the combined weight of the fixture, the aluminum 
mounting plate and the test weight (500, 1000 or 1500 pounds). 

2. Data Collection and Processing 

Three piezoelectric-type accelerometers were utilized during Phase II testing. 
The data collection system captured two seconds of data triggered after MWSM hammer 
release. Through post-processing, each data time record was standardized to one second 
in length triggered off of the initial impulse peak to which an identical negative delay was 
applied for all hammer drops. 

The three accelerometers utilized for Phase II testing were positioned to 
record the acceleration at three specific locations. The accelerometer corresponding to 
Channel 1, Endevco Shock Accelerometer Model 7270-60K, was mounted atop the anvil 
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table. The accelerometer for Channel 3, Model 7270-20K, was mounted atop the fixture’s 
intermediate structure near the center. The accelerometer for Channel 5, Model 7270- 
20K, mounted atop the aluminum plate next to the test weight, provided the base 
excitation input required for all shock spectra calculations. A block diagram of the shock 
data acquisition and processing instrumentation is provided in Figure 8. 

Periodic fixture bolt retightening and visual checks for possible plastic yielding 
or cracking of fixture components were performed throughout testing. A four foot level 
was also used to verify that upper and lower spring beams remained straight, particularly 
after blows from higher hammer heights. The accelerometer for Channel 3 on the fixture’s 
intermediate structure appeared to saturate after the first blow of Test 3 (with a test 
weight of 1500 pounds and a hammer height of 4.5 feet), and was replaced with a spare 
Model 7270-60K accelerometer. Phase II testing was completed with the fixture intact in 


its original condition. 


B. PHASE II DATA REDUCTION 


Of the 180 data files collected for the three accelerometer channels during 60 blows 
of Phase II testing, five data files were not usable due to signal analyzer input voltage 
overload. Of the three accelerometer channels recorded, Channel 5, which measured the 
acceleration at the base of the test weight, provided the “output” of the 2DOF-TDS and, 
therefore, the “input” to the test weight atop the 2DOF-TDS. The test weight, 


representing a combat systems component requiring shock testing, was modeled as an 
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undamped single DOF system subjected to a base excitation for all subsequent absolute 
acceleration shock spectra calculations. 
1. Absolute Acceleration Shock Spectra Determination Using MATLAB” 
The MATLAB” program, originally developed by Cox to calculate shock 
spectra, is contained in Appendix C.1. The absolute acceleration shock spectra were 
obtained through the simulation of a linear time-varying state-space model characterized 


by two equations: 


X()=A(Dx(Y+ Bul) Q) 


IY = C(x) + D(Qu(y) (2) 


where u(t) is the input vector, y(t) is the output vector, and x(t) is the state vector. The 
matrices A(t), Bft), C(t), and D(t) are the state matrix, input matrix, output matrix, and 
direct transmission matrix, respectively (Ogata, 1992, p. 289). The equation of motion for 
an undamped single-DOF system subjected to base excitation was chosen to represent the 


test weight atop the 2DOF-TDS on the MWSM. 


FO +O yf) = -2() (3) 


In Equation (3), y(t) and jt) represent the relative displacement and the 
relative acceleration, respectively, between the top of the 2DOF-TDS and the center of 
mass of the test weight. The single DOF system natural frequency in radians per second is 


represented as @,. The absolute acceleration input, z(t), is the acceleration time history 
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recorded by accelerometer Channel 5 located atop the aluminum plate near the base of the 
test weight. 
In order to write Equation (3) in the general form of the state-space model 


defined in Equations (1) and (2), let: 


ly) | 
x(t) =| (4a) 
and 
TO) 
#8) =| | ao 
yf) = [50] (Sa) 
and 
u(y) =(z(9] (5b) 


and substitute Equations (4) and (5) into Equations (1) and (2), to obtain the state-space 
matrix representation of the undamped single DOF system atop the 2DOF-TDS in 


Equations (6) and (7). 


(sa! _[ o ee fol. 
ice Lye! Leo: olla |*|-z|2! 6) 
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yo 


o An -1[z(1)] (7) 


y(0)=[y(0]=[-0F ol] 


The “Isim” command in the MATLAB® Control System Toolbox was then 
used to calculate the acceleration time response of the undamped single-DOF system to 
input u(t), which contains the absolute acceleration data from accelerometer Channel 5. 


The MATLAB’ code statement: 
[yspec]=Isim(E, F, GG, H, accel, t) (8) 


utilizes the input variables (E, F, GG, H, accel, t) to calculate an output, yspec, which is 
the relative acceleration of the system for a specified system natural frequency calculated 
at each time for which the base acceleration was recorded on Channel 5. The variables, E, 
F, GG and H, correspond to the coefficients of state-space matrix Equations (6) and (7) 
which in turn correspond to the coefficients of Equations (1) and (2), or A, B, C and D, 
respectively. The accel input is the actual recorded accelerometer Channel 5 data 
converted into an ASCII code file. This file contained 2048 consecutive, equally-spaced 
samples of the input acceleration over a sample duration of one second. The t input is the 
one second sample duration divided into discrete sample intervals of 0.4883 milliseconds. 
The relative acceleration yspec was then added to the absolute acceleration 


input accel recorded at the base of the test weight to obtain the absolute acceleration of 
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the single-DOF system, xspec. Since the absolute maximum response of the undamped 


single-DOF system was required, the MATLAB” program code: 


maxspec(i)=max(abs(xspec)) (9) 


stores the largest value of absolute acceleration computed at each system natural 
frequency in the array Maxspec. The absolute acceleration shock spectra for each blow 
was then plotted versus system natural frequency. The frequency band of interest in the 
characterization of the 2DOF-TDS was 0 to 150 Hz stepped in 1 Hz increments. 


2. Resultant Shock Spectra from Phase II Tests 


An illustration of the acceleration time history recorded on Channel 5 for Test 
4 Blow 7 over the resultant shock spectra is presented in Figure 9. Additionally, Figures 
10 through 15 contain the shock spectral plots for the respective six tests. In all tests with 
the exception of Test 1 and Test 5, shock spectra from at least one of the ten blows were 
not plotted due either to abnormal variance in peak amplitude from the norm or to signal 
analyzer input voltage overload during that particular blow. 

While Phase II test results in the form of shock spectra do characterize the 
2DOF Tuned Deck Simulator over six sets of test conditions, these results do not reveal 
the potential advantages of using the 2DOF-TDS in place of the standard fixture. In order 
to appreciate the value of the 2DOF-TDS to shock qualification, a comparison between 
shock spectra generated by the MWSM and standard fixture and those generated in Test 


Phase II by the MWSM and 2DOF-TDS under similar test conditions follows. 
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C. STANDARD FIXTURE SPECTRA VERSUS 2DOF 
TUNED DECK SIMULATOR SPECTRA 


This comparison utilized the acceleration data obtained by the Naval Undersea 
Warfare Center (NUWC), New London, Conn. in a June 1993 test series to characterize 
the MWSM with standard fixture and the acceleration data from Phase II of the 2DOF- 
TDS characterization shock tests. The shock spectra generated by the fixtures were 
compared for test weights of 500 and 1000 pounds. In both cases, the 300 pound 
aluminum plate was attached to the top of each fixture to support the test weight. 


1. Comparison of Shock Spectra with Test Weight of 500 Pounds 


For the 500 pound test, the 2DOF-TDS was configured with three upper 
spring beams and four lower spring beams. Since the 2DOF-TDS weighed more than the 
standard fixture, compliance with Table 1 of MIL-S-901D required that the hammer 
height for the 2DOF-TDS test be one foot higher than that of the standard fixture test. 
The anvil table travel limit was set at three inches for both tests. The accelerometer used 
for calculating shock spectra was mounted on the aluminum plate adjacent to the test 
weight. 

Figure 16 depicts the differences in the respective shock spectra. The first 
spectral peak for the MWSM and standard fixture occurred at 78 Hz with an acceleration 
of 307 g’s. The first peak for the MWSM and 2DOF-TDS occurred at 64 Hz with an 
acceleration of 2017 g’s. Similarly, the standard fixture’s second peak occurred at 162 Hz 
with an acceleration of 180 g’s; the 2DOF-TDS second peak occurred at 102 Hz with an 


acceleration of 1083 g’s. 
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2. Comparison of Shock Spectra with Test Weight of 1000 Pounds 
For the 1000 pound test weight, the center spring beam in the upper tier of the 
2DOF-TDS was removed. The hammer heights for the standard fixture and 2DOF-TDS 
were 2.25 feet and 3.5 feet, respectively. The anvil table travel limit for both tests was 1.5 
inches. 
Figure 17 depicts the differences in the respective shock spectra) The MWSM 
and standard fixture generated shock spectral peaks at 63 Hz (245 g’s), 118 Hz (340 g’s) 
and 184 Hz (772 g’s). The MWSM and 2DOF-TDS spectral peaks occurred at 51 Hz 
(1472 g’s), 90 Hz (551.8 g’s) and 122 Hz (440 g’s). 
3. Summary of Observations 
In both shock spectra comparisons, the 2DOF-TDS tests required higher 
hammer drop heights because the 2DOF-TDS weighed more than the standard fixture. In 
both the 500 and 1000 pound tests, the 2DOF-TDS transmitted shock spectra contained 
initial and secondary peaks of greater amplitude and at lower frequencies than those 
transmitted by the standard fixture. A summary of comparison results for the first and 
second spectral peaks follows: 
«500 POUND TEST WEIGHT: 


=> Amplitude of first spectral peak of 2DOF-TDS was over 6.5 times 


greater than that of standard fixture. 


=> Frequency of first spectral peak of 2DOF-TDS was 14 Hz lower than 


that of standard fixture. 
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=> Amplitude of second spectral peak of 2DOF-TDS was six times greater 
than that of standard fixture. 
=> Frequency of second spectral peak of 2DOF-TDS was 60 Hz lower than 
that of standard fixture. 
¢ 1000 POUND TEST WEIGHT: 
=> Amplitude of first spectral peak of 2DOF-TDS was six times greater 
than that of standard fixture. 
=> Frequency of first spectral peak of 2DOF-TDS was 12 Hz lower than 
that of standard fixture. 
=> Amplitude of second spectral peak of 2DOF-TDS was over 1.5 times 
greater than that of standard fixture. 
=> Frequency of second spectral peak of 2DOF-TDS was 28 Hz lower than 
that of standard fixture. 
It should be noted that the amplitude of the standard fixture’s third spectral peak was over 
1,7 times greater than that of the 2DOF-TDS. The frequency of the third spectral peak of 
the standard fixture was 62 Hz higher than that of the 2DOF-TDS. 
4. Discussion of Results 
According to Clements, increasing the height of the hammer drop raises the 
level of the shock spectrum without changing its shape, while increasing the load weight 
reduces shock severity for frequencies below 200 Hz (Clements, 1972, pp. 94-95). 
Clements also regarded anvil table travel limits as having no consistent effect on the 


generated shock spectra. This limited comparison between shock spectra generated by the 
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standard fixture and the 2DOF-TDS neither supports nor disputes Clements’ description 
of MWSM shock phenomena. The removal of the center spring beam in the 2DOF-TDS 
upper tier and the structural interaction between tiers were not available for Clements to 
characterize. Nonetheless, this comparison does address the issue of which MWSM 
addition, the standard fixture or the 2DOF Tuned Deck Simulator, more realistically 
simulates the actual shock phenomena experienced by deck-mounted equipment. 

If the ultimate objective is to enable the MWSM to reproduce shock spectra 
with at least two dominant spectral peaks at the relatively lower frequencies typically 
transmitted through a ship’s superstructure during shock trials, the 2DOF Tuned Deck 
Simulator approaches this objective more closely than the standard fixture. Incorporating 


the 2DOF-TDS into current shock qualification procedures is the focus of Chapter IV. 
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IV. SHOCK QUALIFICATION AND THE 2DOF 
TUNED DECK SIMULATOR 


The 2DOF Tuned Deck Simulator for the MWSM was designed and built to excite 
combat systems equipment with most of the energy concentrated at two dominant 
resonant frequencies pre-selected from actual ship class shock trial spectral analysis or 
from modal analysis at the equipment’s designated location aboard ship. Cox based his 
2DOF-TDS design on the ability to generate dominant acceleration shock spectral peaks 
at 30 Hz and 80 Hz, respectively (Cox, 1993, pp. 28-29). The lowest frequency 
associated with an initial spectral peak generated in Phase II characterization testing was 
43 Hz (achieved during Test 4 blows). 

This chapter presents the actual frequencies of shock spectral peaks achieved with 
the 2DOF-TDS on the MWSM during Phase II testing and compares these frequencies 
with those obtained using the 2DOF Tuned Deck Simulator’s MATLAB’ and 
MATHCAD* models, respectively. The 2DOF-TDS compatibility with current 


mediumweight shock qualification procedures and its potential integration into MIL-S- 


901D are also explored. 


A. PHASE II ACCELERATION SHOCK SPECTRAL 
PEAKS 


The absolute acceleration shock spectra for blows in Test 1 through Test 6, depicted 


in Figures 10 through 15, reveal at least two (in most cases, three) distinct spectral peaks. 
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The average absolute acceleration shock spectra for each Phase II test and the 
corresponding “modeled” average absolute acceleration shock spectra using a MATLAB” 


model of the 2DOF Tuned Deck Simulator are depicted in Figures 18 through 23. 


1. ‘Effects of Averaging Absolute Acceleration Shock Spectra 

The primary objective of Test Phase II was to increase the accuracy of the 
acceleration shock spectra that best characterized the 2DOF-TDS on the MWSM. By 
increasing the sample size of each test from two hammer drops (in Phase I) to ten drops 
(in Phase II) and averaging the resultant shock spectra for each test, one would expect a 
reduction in the average random noise in the data and an increase in spectral peak 
definition. 

In almost all tests, however, the initial spectral peaks survived the averaging 
process intact, but the second spectral peaks were either greatly diminished or vanished 
completely after averaging. For example, Figure 10 depicts the absolute acceleration 
shock spectra for the ten drops in Test 1. A second spectral peak at 102 Hz, clearly 
evident in Figure 10, is conspicuously missing from the average shock spectra plot in 
Figure 18. Since all Phase II acceleration data files were standardized in post-processing 
to a one second time record triggered off of the initial impulse peak and adjusted with the 
same negative time delay, trigger initiation during data collection was not a contributing 
factor in the reduction or loss of second spectral peaks upon averaging. Instead, slight 
phase shifts are evident in Figures 10 through 15 between the second spectral peaks of 
various drops in each test. When averaging the absolute acceleration values for the 


respective spectra in each test, any “out-of-phase” second spectral peaks could effect a 
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reduction in second peak amplitude. Possible causes of the phase shifts include: non- 
uniform bolt tightening on the 2DOF-TDS between hammer drops; and slight weight shifts 
of the 500 pound lead weights representing the test object atop the 2DOF-TDS. 

2. “Modeled” Average Absolute Acceleration Shock Spectra 


In addition to displaying the “actual” average absolute acceleration shock 
spectra for each Phase II test, Figures 18 through 23 also present the “modeled” average 
absolute acceleration shock spectra. The MATLAB’ model, used for the comparison 
between “actual” average and “modeled” average, was originally developed by Cox (Cox, 
1993, pp. 109-112). This model was updated with the actual 2DOF-TDS tier weights and 
modified to input actual Phase II acceleration data (Channel 1) obtained at the top of the 
anvil table instead of the half-sine wave vertical acceleration impulse of one millisecond 
duration (used previously by Clements and Corbell). Cox estimated the value of the 
critical damping ratio (zeta) for the 2DOF-TDS to be between three and five percent 
(Cox, 1993, p.23). The MATLAB® model was again updated after a logarithmic 
decrement calculation for zeta using acceleration time histories obtained during Phase | 
testing yielded an average zeta of 0.035. 

The MATLAB?’ program that modeled 2DOF Tuned Deck Simulator behavior 
under the same initial conditions as Phase II Test 1 and produced Figure 18 is provided in 
Appendix C.2. The absence of a second spectral peak in the “modeled” average spectra 


for Test 4 (Figure 21) indicated that the “modeled” average spectra were also vulnerable 


to the effects of “out-of-phase” second spectral peaks. 
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3. Determining System Natural Frequencies Using MATHCAD® 
A MATHCAD?® worksheet was utilized by Cox to determine the achievable 
system natural frequencies for various upper tier weights and 2DOF-TDS configurations 
(Cox, 1993, pp. 106-107). The MATHCAD” worksheet for the first and second system 


natural frequencies generated under the conditions of Test 1 is provided in Appendix C.3. 
The 2DOF-TDS in Test 1 was configured with all upper and lower spring 
beams installed and with a 500 pound lead weight atop the aluminum support plate. The 
spring stiffness per spring beam was calculated using the beam bending formula for a 
simply supported beam with two equal concentrated forces symmetrically placed. The 
upper tier stiffness and lower tier stiffness were obtained by multiplying the spring stiffness 
per beam by the respective number of spring beams per tier. After entering the actual 
weights of the upper tier, intermediate structure and anvil table, the tier natural 


frequencies and system natural frequencies were computed. 


B. FREQUENCIES OF SHOCK SPECTRAL PEAKS 


1. Methods Utilized 


The frequencies of the absolute acceleration shock spectral peaks in Test 
Phase II were determined using three methods, specifically: by graphically estimating 
spectral peak locations in Figures 10 through 15; by averaging the respective shock 
spectra of each test and determining spectral peak locations from the average shock 
spectra plots in Figures 18 through 23; and by averaging the MATLAB” “modeled” shock 


spectra of each test and determining spectral peak locations from the “modeled” shock 
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spectra plots also in Figures 18 through 23. Comparative summaries of spectral peak 
locations, including the system natural frequencies from MATHCAD® worksheets, for 
Tests 1 through 6 are provided in Tables 2 through 13. 

Due to the effects of “out-of-phase” second spectral peaks on shock spectra 
averaging, the most accurate spectral peak locations were obtained via the graphical 
estimation method using Figures 10 through 15. The spectral peak locations obtained 
from the “modeled” average shock spectra for each test were included in the summary 
tables to illustrate how closely the MATLAB® model’s average spectral peak locations 
matched those obtained via graphical estimation. Since the input to the MATLAB” model 
(in Appendix C.2) was acceleration recorded on the anvil table, the “modeled” shock 
spectra were considered less accurate than the actual shock spectra based on acceleration 
recorded at the base of the test weight atop the 2DOF-TDS. 

The two system natural frequencies obtained via the MATHCAD* worksheet 
(in Appendix C.3) were also included in Tables 2 through 13. These frequencies 
characterized the undamped free vibration of the system from the hammer impact point 
on the anvil table to the center of mass of the test weight on the 2DOF-TDS. Cox 
demonstrated the usefulness of his MATHCAD® worksheet by determining the range of 
system natural frequencies attainable for a specified equipment weight atop the 2DOF 


Tuned Deck Simulator (Cox, 1993, pp. 46-47). 
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2. Comparison Between Methods 


Some interesting trends are apparent among the frequencies of shock spectral 
peaks displayed in Tables 2 through 13. For each of the six tests, the frequencies of the 
first spectral peaks, determined via the four respective methods, were practically identical. 
Although a second spectral peak was apparent in the average shock spectra for four of the 
six tests (depicted in Figures 18 through 23), the amplitude of each second peak was 
significantly lower than that of the corresponding actual shock spectra (in Figures 10 
through 15). This reduction in average acceleration amplitude was attributed to the phase 
differences between various second spectral peaks and their effect on spectra averaging. 

Regarding the spectral peak locations from the graphical estimate method as 
the most accurate, the proximity of second spectral peak locations relative to those 
obtained via the graphical estimate method revealed two other trends. Specifically, the 
second spectral peaks from the MATLAB® “modeled” average method (with the 
exception of Test 4), occurred an average of 9.4 Hz higher than those from the graphical 
estimate method. The MATHCAD”-derived second natural frequency of the 2DOF-TDS 
was an average of 12.67 Hz greater than the frequency of the second spectral peak from 
the graphical estimate method for all six tests. 

Although this characterization of the 2DOF Tuned Deck Simulator and 
validation of its MATLAB® model were limited to six tests on the MWSM, the 2DOF- 
TDS can be integrated into the mediumweight shock qualification specifications of MIL- 


S-901D. 
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C. 2DOF TUNED DECK SIMULATOR 
COMPATIBILITY WITH MIL-S-901D 


1. __ Basic Requirements 

In order to be compatible with the current test schedule and fixture 
requirements for mediumweight shock testing in MIL-S-901D, the total weight on the 
anvil table must not exceed 7400 pounds; and the 2DOF-TDS must accommodate the 
standard inclined fixture in order to simulate shipboard shock loadings in the athwartship 
direction (MIL-S-901D, 1989, pp. 16-17). The 2DOF-TDS with all spring beams 
installed (and the aluminum plate removed) weighs 2675 pounds. The standard inclined 
fixture, which can be easily mounted onto the 2DOF-TDS at the two I-beams which cross 
the upper tier spring beams, weighs 715 pounds. Since three out of the six hammer blows 
must be conducted with the tested item in the inclined position, the test item’s weight is 
limited to 4010 pounds when mounted on the 2DOF-TDS. 

MIL-S-901D further stipulates that fixtures intended to simulate deck- 
mounted conditions shall possess a fundamental response frequency (with test item 
mounted) not lower than 25 Hz and that any plastic yielding or cracking of fixtures is 
unacceptable (MIL-S-901D, 1989, p. 13). The lower limit on 2DOF-TDS fundamental 
response frequency is not considered restrictive as the 2DOF-TDS design’s target 
frequency was 30 Hz, and the lowest frequency attained during both Test Phases I and II 


was 43 Hz. 
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2. Maximum Bending Stress Versus Yield Strength 
Regarding plastic yielding or cracking of the 2DOF-TDS, the spring beams 
were made of carbon steel with a yield strength of 50 ksi. All other structural elements (in 
Figure 6) were made of mild steel with a yield strength of 36 ksi. A total of 84 zinc-plated 
steel socket head bolts (1-8UNC x 4.25 inches) were used in the construction of the 
2DOF-TDS 
Cox predicted a maximum bending stress in the spring beams of 47 ksi based 


on his MATLAB’ model of the 2DOF-TDS, loaded to its full weight capacity and excited 


with Clements’ maximum expected anvil acceleration of 200 g’s for a five foot hammer 
drop (Cox, 1993, pp. 39-43). During Test Phase I, however, after a hammer drop of 3.75 
feet with a test weight of 1500 pounds atop the 2DOF-TDS (in “2/4” configuration), a 
deflection of approximately '/; of an inch at the center of the upper spring beam nearest 
the MWSM hammer was observed. The effected spring beam was subsequently 
straightened, and testing continued without incident. 

Although this spring beam deflection was attributed to a slight lateral shift 
toward the hammer end of the MWSM of the test weights (three 500 pound lead cylinders 
stacked vertically on the aluminum plate) prior to the 3.75 foot hammer drop, the total 
weight on the anvil table for this Phase I test was 4499 pounds, equating to 60 percent of 
anvil table capacity. The 3.75 foot height for this test equates to 68 percent of the 
maximum hammer height prescribed by MIL-S-901D. No beam deflection or test weight 


lateral movement was observed during Test Phase II; however, the hammer height for the 
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ten drops in Phase II Test 4 (with a test weight of 1500 pounds atop the 2DOF-TDS in 
“2/4” configuration), was only 2.25 feet. 

Since the possibility that at least one of the 2DOF-TDS spring beams could 
plastically deform during a shock test at greater than 70 percent of MWSM capacity has 
not been ruled out, a half-bridge strain gage configuration (with two active arms for the 
respective I-beam flange surfaces and with temperature compensation) could be mounted 
at the center of both an upper tier and a lower tier spring beam. The 2DOF-TDS in the 
“2/4” and “3/4” configurations could then be incrementally loaded to capacity and shock 
tested at MIL-S-901D specified hammer heights in order to determine the maximum 
bending stress in the spring beams. Should the maximum bending stress exceed the 50 ksi 
yield strength of carbon steel and/or plastic deformation occur during shock testing, some 
restrictions on equipment test weight and hammer height when using the 2DOF-TDS may 
be required. To avoid these restrictions, the spring beams could be made from steels with 
higher yield strengths. 

3. Compatibility Achieved 

Based purely on the results of Test Phase Il (i.e., barring the possibility of 
spring beam deformation under other than Phase II test conditions), the 2DOF Tuned 
Deck Simulator meets the current MWSM test schedule and fixture requirements of MIL- 
$-901D. 

For the 2DOF Tuned Deck Simulator to be fully compatible with MIL-S- 
901D and utilized to its full potential, the mediumweight shock test criteria for deck- 


mounted equipment should reflect the ultimate design goal for the 2DOF-TDS; that is to 
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excite combat systems equipment with most of the energy concentrated at two dominant 
resonant frequencies pre-selected from actual ship class shock trial spectral analysis or 
from modal analysis at the equipment’s designated location aboard ship. Only under this 
context would the 2DOF-TDS offer significant advantages over the standard fixture in 
mediumweight shock qualification. A proposal for the integration of the 2DOF-TDS into 


MIL-S-901D is contained in Chapter V. 


35 


Vv. CONCLUSION 


This experimental investigation into the dynamic response of the 2DOF Tuned Deck 
Simulator has proved that the 2DOF-TDS can be successfully integrated into the U.S. 
Navy’s mediumweight shock qualification program. Unlike the standard fixture for the 
MWSM, the 2DOF-TDS demonstrated the potential to excite combat systems equipment 
with most of the energy concentrated at two dominant resonant frequencies pre-selected 
from actual ship class shock trial spectral analysis or from modal analysis at the 


equipment’s designated location aboard ship. 


A. PROJECTED UTILIZATION OF THE 2DOF 
TUNED DECK SIMULATOR 


With the emphasis now on transmitting a high acceleration impulse at both the first 
and second resonant frequencies (characteristic of a specific shipboard location) into the 
base of a combat systems component via the 2DOF Tuned Deck Simulator, the current 
test schedule for the MWSM in Table 1 of MIL-S-901D can be simplified (MIL-S-901D, 
1989, p. 17). Instead of requiring three groups of two blows each at various hammer 
heights and anvil table travel limits, the test schedule could require just two groups of two 
blows each, in which the hammer height and anvil table travel limit were based only on the 
desired primary and secondary resonant frequencies and on the total weight on the anvil 


table. 


36 


For both Group I blows, the test item would be mounted in its normal attitude atop 
the 2DOF-TDS to simulate shock loadings in the vertical direction. For both Group II 
blows, the test item would be mounted on the standard inclined fixture atop the 2DOF- 
TDS to simulate shock loadings in the athwartship direction. The test item would be 
operated in its normal mode during the first blow and in its alternate mode during the 


second blow of each group. 


B. PROPOSED ADDENDUM TO MIL-S-901D 


1. Proposed Test Schedules 


A proposed addendum to the mediumweight shock qualification procedures in 
MIL-S-901D based solely on Phase II test data is provided in Tables 14 and 15. The test 
schedule in Table 14 applies when the initial and second spectral peaks occur within the 
ranges of 45 - 55 Hz and 90 - 100 Hz, respectively. The test schedule in Table 15 applies 
when the initial and second spectral peaks occur within the ranges of 55 - 65 Hz and 90 - 
100 Hz, respectively. Directions for computing the “Total Weight on the Anvil Table” are 
listed after Table 15. 


2. Need for Expanded Coverage 
Before the proposed addendum could be integrated into MIL-S-901D as a test 
schedule for the MWSM, Tables 14 and 15 would have to be expanded to cover greater 
ranges of equipment weights and spectral peak frequencies. More diverse characterization 
testing of the 2DOF-TDS using multiple “tuning” features under a wider range of test 


weights would be necessary. 
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Since the MATHCAD* worksheet (in Appendix C.3) enables the user to 
“tune” the 2DOF-TDS model by removing the center spring beam from the upper tier, by 
altering the load application points on the upper spring beams, or by adding ballast to the 
intermediate structure, the MATHCAD® worksheet would be useful in the formulation of 
a prediction-based test plan for a final phase of 2DOF-TDS characterization testing. 

3. Environmental Issues 

The future integration of the 2DOF-TDS into mediumweight shock 
qualification procedures may be driven more by environmentalists who oppose ship shock 
trials than by NAVSEA proponents of more accurate shock qualification methods. 
According to a 14 February 1994 article in the San Francisco Chronicle, an ad hoc group 
of environmentalists are exploring legal challenges to the shock trial of DDG-53, originally 
scheduled to occur in March 1994 in an offshore test range whose border was 
approximately six miles from the Channel Islands National Marine Sanctuary off the 
California coast (Martin, 1994, p.13). 

The U.S. Navy’s commitment to conducting surface ship shock trials on one 
ship of each shock hardened class remains steadfast (NAVSEA , 1989, pp.2-4). The 
Commanding Officer of an AEGIS guided missile cruiser that underwent shock trials in 
1984 commented that ship shock trials are cost effective and minimize future risk and 
casualties (Anderson, 1987, p. 39). 

If the courts side with the environmentalists, and the Navy is unable to 
conduct shock trials in areas deemed safe by environmental groups, more accurate shock 


qualification may become the only feasible alternative for the near term. Since equipment 
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operability following shock is verified after mediumweight shock test blows, the 2DOF 
Tuned Deck Simulator and MWSM could provide the more accurate and comprehensive 


method of mediumweight shock qualification. 


C. RECOMMENDATIONS 


While the 2DOF-TDS is compatible with the mediumweight shock qualification 
requirements and procedures of MIL-S-901D, an additional test phase on the MWSM is 
necessary to completely characterize the 2DOF-TDS under test weights and “tuning” 
configurations other than those previously tested. The MATHCAD’ worksheet in 
Appendix C.3 would be useful in formulating a prediction-based test schedule for future 
characterization testing. 

As part of a future test phase, strain gages should be mounted on upper and lower 
tier spring beams of the 2DOF Tuned Deck Simulator to determine the maximum bending 
stress in the spring beams during shock testing. Should the maximum bending stress 
exceed the 50 ksi yield strength of carbon steel and/or plastic deformation occur during 
shock testing, some restrictions on equipment test weight and hammer height when using 
the 2DOF-TDS may be required. To avoid these restrictions, the spring beams could be 
made from steels with higher yield strengths. 

In addition to expanding the characterization of the 2DOF-TDS to complete the 
proposed addendum to MIL-S-901D, one possible direction for new research is the 
development of a MWSM fixture capable of generating a first spectral peak between 25 


and 40 Hz. 
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Once completed and approved, the addendum to MIL-S-901D for the utilization of 
the 2DOF Tuned Deck Simulator should provide a quantum improvement in the 
simulation of a ship’s actual response to an underwater explosion for the pre-acceptance 
shock qualification of combat systems equipment. This improvement promises to enhance 


the capacity of a warship to absorb damage and still maintain its mission integrity. 
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FIGURES 














APPENDIX A 





Figure 1. The Navy High-Impact Shock Machine for Mediumweight (MWSM). 


The Dotted Line Shows Hammer Path (Courtesy of Clements). 
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DDG-51 Pretrial Prediction Versus MWSM with Standard Fixture 


NOTES: (1) Prediction for SPY-1D Beam Programmer (1000 Ibs) 
(2) Weight on MWSM Standard Fixture = Approx. 1300 Ibs 
(3) MWSM Hammer Height = 2.25 ft; Anvil Travel Limit = 3 in 
(4) 400 Hz Bessel Filter Used in MWSM Test 


Acceleration Time History for SPY-1D Beam Programmer 
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Figure 2. DDG-51 Pretrial Predicted Acceleration Time History Versus MWSM 
with Standard Fixture Actual Acceleration Time History. 
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DDG-51 Pretrial Prediction Versus MWSM with Standard Fixture 


NOTE: (1) Prediction for SPY-1D Beam Programmer (1000 Ibs) 
(2) Weight on MWSM Standard Fixture = Approx. 1300 Ibs 
(3) MWSM Hammer Height = 2.25 ft; Anvil Travel Limit = 3 in 
(4) 400 Hz Bessel Filter Used in MWSM Test 
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Figure 3. DDG-51 Pretrial Predicted Absolute Acceleration Shock Spectra Versus 
MWSM with Standard Fixture Actual Absolute Acceleration Shock Spectra. 
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M1 ee as 
M, = Upper Tier Mass 
Includes: Equipment Mass 


Mounting Hardware 
C1 «1 Effective Spring Mass 


M, = Lower Tier Mass 


a x2 Includes: Intermediate. Mass 
M2 Effective Spring Mass 


K = Tier Stiffness 


C2 C = Tier Structural Damping 


Figure 4. Model of 2DOF Tuned Deck Simulator as Two Mass-Spring-Damper 
Systems Coupled Together (Courtesy of Cox). 
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Figure 5. 1/4 Scale Model of the 2DOF Tuned Deck Simulator for the MWSM. 
(Courtesy of Cox). 
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Figure 6. Bill of Mater 


Silver Springs, Md.). 


Drawing (Courtesy of Naval Surface Warfare Center, 
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Upper Tier Aluminum Test Weight Equipment Specific 


Spring Beams (3) Plate (500 Ibs) Mounting Beams (2) 
“s |, 
ye :, 






Intermediate Lower Tier Anvil Table with Existing 
Structure Spring Beams (4) Shipbuilding Rails (2) 


Figure 7. Photograph of 2DOF Tuned Deck Simulator with Aluminum Mounting 
Plate and 500 Pound Test Load atop the MWSM Anvil Table. 
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; ENDEVCO SHOCK HEWLETT-PACKARD 
; ACCELEROMETERS SIGNAL ANALYZER 


Channel Model Model 35665 


1 7270-60K 
2 7270-60K 
3 7270-20K 





PRISM 


Portable Personal Computer 
(486/33 Mhz) 






ENDEVCO SIGNAL 
CONDITIONER 
Model 4423 


HEWLETT-PACKARD 
Printer (Laserjet 4) 


Figure 8. Block Diagram of Test Phase II Shock Data Acquisition and Processing 
Instrumentation. 
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MWSM 2DOF Tuned Deck Simulator - Test 4 Blow 7 


Spring Beam Configuration: 2 Upper / 4 Lower 
Test Weight on Deck Simulator: 1500 Ibs 

Total Weight on Anvil Table: 4389 Ibs 

MWSM Hammer Ht: 2.25 ft (Anvil Travel = 3 in) 


Acceleration Time History (Channel 5) 
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Figure 9. Acceleration Time History (Channel 5) and Absolute Acceleration Shock 
Spectra for Test 4 Blow 7. 
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ABSOLUTE ACCELERATION SHOCK SPECTRA FOR TEST 1 (10 DROPS) 
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Figure 10. Absolute Acceleration Shock Spectra for Test 1 (10 Drops). 
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ABSOLUTE ACCELERATION SHOCK SPECTRA FOR TEST 2 (8 DROPS) 
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Figure 11. Absolute Acceleration Shock Spectra for Test 2 (8 Drops). 


31 


ABSOLUTE ACCELERATION SHOCK SPECTRA FOR TEST 3 (9 DROPS) 
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Figure 12. Absolute Acceleration Shock Spectra for Test 3 (9 Drops). 
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ABSOLUTE ACCELERATION SHOCK SPECTRA FOR TEST 4 (9 DROPS) 
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Figure 13. Absolute Acceleration Shock Spectra for Test 4 (9 Drops). 
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ABSOLUTE ACCELERATION SHOCK SPECTRA FOR TEST 5 (10 DROPS) 
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Figure 14. Absolute Acceleration Shock Spectra for Test 5 (10 Drops). 
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ABSOLUTE ACCELERATION SHOCK SPECTRA FOR TEST 6 (9 DROPS) 
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Figure 15. Absolute Acceleration Shock Spectra for Test 6 (9 Drops). 


Standard Fixture Versus 2DOF Deck Simulator (3/4 Configuration) 


NOTE: (1) Aluminum Plate (300 ibs) Supports Test Weight. 
(2) Test Weight atop Aluminum Plate = 500 Ibs. 
(3) Hammer Ht (HH) and Anvil Travel (AT) per MIL-S-901D. 
(4) 400 Hz Bessel Filter Used in Standard Fixture Test. 
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Figure 16. Standard Fixture Versus 2DOF Tuned Deck Simulator (3/4 
Configuration). 
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Standard Fixture Versus 2DOF Deck Simulator (2/4 Configuration) 


NOTE: (1) Aluminum Plate (300 Ibs) Supports Test Weight. 
(2) Test Weight atop Aluminum Plate = 1000 Ibs. 
(3) Hammer Ht (HH) and Anvil Travel (AT) per MIL-S-901D. 
(4) 400 Hz Bessel Filter Used in Standard Fixture Test. 


Shock Spectra for Standard Fixture (HH=2.25 ft; AT=1.5 in) 
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Figure 17, Standard Fixture Versus 2DOF Tuned Deck Simulator (2/4 
Configuration). 
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2DOF Tuned Deck Simulator - Test 1 

Spring Beam Configuration: 3 Upper / 4 Lower 
Test Weight on Deck Simulator: 500 Ibs 

Total Weight on Anvil Table: 3494 Ibs 

MWSM Hammer Et: 3 ft (Anvil Travel = 3 in) 
10 of 10 Hammer Drops Used for Average. 


Average Absolute Acceleration Shock Spectra 
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Figure 18. Average Absolute Acceleration Shock Spectra Versus “Modeled” 
Average Acceleration Shock Spectra for Test 1. 
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2DOF Tuned Deck Simulator - Test 2 

Spring Beam Configuration: 3 Upper / 4 Lower 
Test Weight on Deck Simulator: 1000 Ibs 

Total Weight on Anvil Table: 3994 Ibs 

MWSM Hammer Et: 3.5 ft (Anvil Travel = 1.5 in) 
8 of 10 Hammer Drops Used for Average. 
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Figure 19. Average Absolute Acceleration Shock Spectra Versus “Modeled” 
Average Acceleration Shock Spectra for Test 2. 
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2DOF Tuned Deck Simulator - Test 3 

Spring Beam Configuration: 3 Upper / 4 Lower 
Test Weight on Deck Simulator: 1500 Ibs 

Total Weight on Anvil Table: 4494 Ibs 

MWSM Hammer Ht: 4.5 ft (Anvil Travel = 1.5 in) 
9 of 10 Hammer Drops Used for Average. 
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Figure 20. Average Absolute Acceleration Shock Spectra Versus “Modeled” 
Average Acceleration Shock Spectra for Test 3. 
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2DOF Tuned Deck Simulator - Test 4 

Spring Beam Configuration: 2 Upper / 4 Lower 
Test Weight on Deck Simulator: 1500 Ibs 

Total Weight on Anvil Table: 4389 ibs 

MWSM Hammer Ht: 2.25 ft (Anvil Travel = 3 in) 
9 of 10 Hammer Drops Used for Average. 
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Figure 21. Average Absolute Acceleration Shock Spectra Versus “Modeled” 
Average Acceleration Shock Spectra for Test 4. 
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2DOF Tuned Deck Simulator - Test 5 

Spring Beam Configuration: 2 Upper / 4 Lower 
Test Weight on Deck Simulator: 1000 Ibs 

Total Weight on Anvil Table: 3889 lbs 

MWSM Hammer Ht: 3.5 ft (Anvil Travel = 1.5 in) 
10 of 10 Hammer Drops Used for Average. 
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Figure 22, Average Absolute Acceleration Shock Spectra Versus “Modeled” 
Average Acceleration Shock Spectra for Test 5. 
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Ave. Accel. 
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2DOF Tuned Deck Simulator - Test 6 

Spring Beam Configuration: 2 Upper / 4 Lower 
Test Weight on Deck Simulator: 500 Ibs 

Total Weight on Anvil Table: 3389 Ibs 

MWSM Hammer Ht: 1.75 ft (Anvil Travel = 3 in) 
9 of 10 Hammer Drops Used for Average. 
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Figure 23. Average Absolute Acceleration Shock Spectra Versus “Modeled” 
Average Acceleration Shock Spectra for Test 6. 
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APPENDIX B: TABLES 


Appendix B.1: Phase II Test Schedule 


TABLE 1. PHASE Il TEST SCHEDULE FOR 2DOF-TDS ON MWSM 






Test Test Weight ' Spring Beam | Height of _ Anvil Table | 
Number on 2DOF- : Configuration : Hammer — Travel (in) 


TDS _ Drop i 
1 500 3 upper/4 lower 3 3 
2 1000 3/4 : 3.5 : 1.5 i 
3 1500 3/4 4.5 1.5 
Remove Center Spring Beam From Upper Tier | 
4 1500 2/4 / 2.25 : 3 i 
5 1000 2/4 a8 1.5 
6 500 2/4 1.75 i 3 








Appendix B.2: Results of Test #1 


TABLE 2. TEST 1: INITIAL CONDITIONS 


Weight on 2DOF Tuned Deck Simulator 500 lbs 
Total Weight on Anvil Table 3494 Ibs 


Spring Beam Configuration 3 upper / 4 lower 
MWSM Hammer Height 3 ft 
Anvil Table Travel Limit 3 in 





TABLE 3. TEST 1: ABSOLUTE ACCELERATION SHOCK SPECTRAL PEAK 
FREQUENCIES 


eae 2nd Peak | 3rd Peak 
Graphical Estimate before Averaging 
10 of 10 Drops, 
Actual Average (10 of 10 Drops) |. 64 Az | 
Modeled Average (10 of 10 Drops) 
System Natural Frequencies from Mathcad” 


Worksheet __ 
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Appendix B.3: Results of Test #2 


TABLE 4. TEST 2: INITIAL CONDITIONS 


1000 Ibs 
3994 Ibs 
3 upper / 4 lower 
3.5 ft 
1.5 in 


Weight on 2DOF Tuned Deck Simulator 
Total Weight on Anvil Table 


Spring Beam Configuration 
MWSM Hammer Height 
Anvil Table Travel Limit 





TABLE 5. TEST 2; ABSOLUTE ACCELERATION SHOCK SPECTRAL PEAK 
FREQUENCIES 







3rd Peak | 







IL comune éVodeled oe (8 of 10 Drops. 

| ~ System Natural Frequencies from Mathcad® _ 
q 

__Worksheet _ 
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Appendix B.4: Results of Test #3 


TABLE 6. TEST 3: INITIAL CONDITIONS 


Weight on 2DOF Tuned Deck Simulator 1500 Ibs 
Total Weight on Anvil Table 4494 lbs 


Spring Beam Configuration 3 upper / 4 lower 
MWSM Hammer Height 4.5 ft 
Anvil Table Travel Limit 1.5 in 





TABLE 7, TEST 3: ABSOLUTE ACCELERATION SHOCK SPECTRAL PEAK 
FREQUENCIES 


| 3rd Peak. 
134 Hz 












System Natural Frequencies from Mathcad* 
Worksheet 
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Appendix B.5: Results of Test #4 


TABLE 8. TEST 4: INITIAL CONDITIONS 


Weight on 2DOF Tuned Deck Simulator 1500 Ibs 
Total Weight on Anvil Table 4389 lbs 


Spring Beam Configuration 2 upper / 4 lower 
MWSM Hammer Height 2.25 ft 
Anvil Table Travel Limit 3 in 





TABLE 9. TEST 4: ABSOLUTE ACCELERATION SHOCK SPECTRAL PEAK 
FREQUENCIES 


3rd Peak 








Graphical Estimate before Averaging 
_(9 of 10 Drops) _ 
~ Actual Average (9 of 10 Drops) 
Modeled Average (9 of 10 Drops) - 
ural Frequencies from Mathcad — 
Worksheet 
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Appendix B.6: Results of Test #5 


TABLE 10. TEST 5: INITIAL CONDITIONS 


Weight on 2DOF Tuned Deck Simulator 1000 Ibs 
Total Weight on Anvil Table 3889 Ibs 


Spring Beam Configuration 2 upper / 4 lower 
MWSM Hammer Height 3.5 ft 
Anvil Table Travel Limit 1.5 in 





TABLE 11. TEST 5: ABSOLUTE ACCELERATION SHOCK SPECTRAL PEAK 
FREQUENCIES 








2nd Peak | 3rd Peak 













_| dst Peak | 





~~ Graphical Estimate before Averaging 


_ _(10 of 10 Drops) : 
Actual Average (10 of 10 Drops) |, SI Hz |) 91 Hz | 122 














System Natural Frequencies from Mathcad® 
i 
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Appendix B.7: Results of Test #6 


TABLE 12, TEST 6: INITIAL CONDITIONS 


Weight on 2DOF Tuned Deck Simulator 500 Ibs 
Total Weight on Anvil Table 3389 Ibs 


Spring Beam Configuration 2 upper / 4 lower 
MWSM Hammer Height 1.75 ft 
Anvil Table Travel Limit 3 in 





TABLE 13. TEST 6: ABSOLUTE ACCELERATION SHOCK SPECTRAL PEAK 
FREQUENCIES 









Ist Peak | 2nd Peak | 3rd Peak 




































Graphical Estimate before Averaging 58 Hz T 
2. of 10 Drops) iene ceed beet Leena teat tet 2 an 
( _ Actual Average (9 of 10 Drops) _ 58 Hz 90 Hz jlllHz | 
sunwnsnendVdodeled Average (9 of 10 Drops) | sok d8HZ. | 99 Hz | INH | 
“System Natural Frequencies from Mathcad® «57.7 Hz | 105.4 Hz - 
Worksheet | 
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Appendix B.8: Proposed Addendum to MIL-S-901D 


TABLE 14, FOR INITIAL SPECTRAL PEAK BETWEEN 45 - 55 HZ AND 
SECOND PEAK BETWEEN 90 -100 HZ 


Total Weight on Spring Beam Height of Hammer | Anvil Table Travel 
Anvil Table Configuration Drop (feet) (inches) 
| (pounds) __—|_—s(upperflower) | 


aaa ae [eT © CE 
oe ETE SE ON: 
7500 










TABLE 15. FOR INITIAL SPECTRAL PEAK BETWEEN 55 - 65 HZ AND 
SECOND PEAK BETWEEN 90 -100 HZ 


















Total Weight on Spring Beam Height of Hammer | Anvil Table Travel 
Anvil Table Configuration Drop (feet) (inches) 


___(pounds) _|_ _(upper/ower) _ per/lower 


(ae ee (ene 1 nr ee ae ee 








3500 
4000 


NOTES: 

1. Group I consists of two blows with the equipment mounted on the 2DOF Tuned 
Deck Simulator. The equipment shall be operated in its normal mode for the first 
blow and in its alternate mode for the second blow. 

2. Group II consists of two blows with the equipment mounted on the Standard 
Inclined Fixture atop the 2DOF-TDS. The equipment shall be operated in its normal 


mode for the first blow and in its alternate mode for the second blow. 


3. With 2DOF Tuned Deck Simulator in “2 / 4” Configuration, the “Total Weight on 
Anvil Table” equals 2570 plus weight of test item in pounds. 


71 


4, With 2DOF Tuned Deck Simulator in “3 / 4” Configuration, the “Total Weight on 
Anvil Table” equals 2675 plus weight of test item in pounds. 


5. Include the weight of the Standard Inclin; Fixture (715 pounds) in the “Total 
Weight on Anvil Table” for all Group II blo: s. 
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APPENDIX C: MATLAB® PROGRAMS AND 
MATHCAD® WORKSHEET 


Appendix C.1: Model APPC1.M 


%THE FOLLOWING MATLAB PROGRAM CALCULATES THE ABSOLUTE ACCELERATION 
%SHOCK SPECTRA FOR TEST 4 BLOW 7 GIVEN ABSOLUTE ACCELERATION DATA 
(FROM ACCELEROMETER CHANNEL 5) OVER A ONE SECOND TIME DURATION 
(CONTAINING 2048 DATA SAMPLES). 
% 
% 
t=linspace(0,1,2048); % Time Interval of 1 second containing 2048 data points. 
load t4b1c5s7.asc 
accel1=accel1(:); 
accel1=t4b1c5s7(1:2048); % Store absolute accel. data in array “accel1.” 
% 
%Calculate the Absolute Acceleration Shock Spectra. 
% 


NF=150; % Number of Frequencies 

DF=1; % Frequency Increment 

SF=1; % Start Frequency 
Freq=[SF:DF:NF]; % Frequency Vector 

wn=2*pi*Freq; % Convert Hz to radians per second. 
F=(0;-1]; % State Space Matrices 

GG=(1,0]; 

H=[-1]; 

for i=1:NF % Step thru natural frequencies. 


E=[0,1;-(wn(i)*2),0]; 
GG=[wn(i)*2,0]; 


[yspect}=Isim(E,F,GG,H,accel1 ,t); % Computes relative accel. 
xspect = yspeci + accel1; % Computes absolute accel. of single 
DOF system. 
maxspec1(i)=max(abs(xspect)); % Stores largest value of absolute 
accel. computed at each natural 
frequency. 
end 
figure(1) 


subplot(31 1), axis(off) 

text(0.1,1,,MWSM 2DOF Tuned Deck Simulator - Test 4 Blow 7') 
text(0.1,0.6,'Spring Beam Configuration: 2 Upper / 4 Lower’) 
text(0.1,0.4,"Test Weight on Deck Simulator: 1500 Ibs’) 
text(0.1,0.2,'Total Weight on Anvil Table: 4389 Ibs’) 
text(0.1,0.0,,MWSM Hammer Ht: 2.25 ft (Anvil Travel = 3 in)') 
subplot(312),plot(t,accel1),grid 

title(Acceleration Time History (Channel 5)’) 

xlabel(‘Time (sec)’) 

ylabei(Accel. (g)’) 
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axis({O 1 -80 80)) 

subplot(313), plot(Freq,maxspec1),grid 
title(Absolute Acceleration Shock Spectra’) 
xlabel(‘Frequency (Hz)} 

ylabel(‘Accel. (g)’) 

axis([0 150 0 800)) 

set(1,‘paperpos’,[1.25 2 6 8]) 

print 

freq=Freq'; 

cd. 
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Appendix C.2: Model APPC2.M 


%THE FOLLOWING PROGRAM CALCULATES THE AVERAGE ABSOLUTE 
%ACCELERATION SHOCK SPECTRA AND THE 2DOF TUNED DECK SIMULATOR 
%MODEL'S AVERAGE ABSOLUTE ACCELERATION SHOCK SPECTRA FOR TEST 
%ONE OF PHASE Il CHARACTERIZATION TESTING. THE ORIGINAL VERSION 
%OF THIS MATLAB CODE WAS PROVIDED BY COX (COX, 1993, PP. 109-112). 


t=linspace(0,1,2048); % TIME INTERVAL OF 1 SECOND CONTAINING 2048 DATA 


POINTS. 
% 
%LOAD=500; % LOAD = TEST LOAD (500, 1000 or 1500 LBS). 
Ynk1=3; % nk1 = NUMBER OF SPRING BEAMS (2 OR 3) IN UPPER TIER. 
Yonk2=4; % nk2 = NUMBER OF SPRING BEAMS IN LOWER TIER. 
cd coxfile 


load tib2c5s1.asc % LOAD DATA FILES FOR 10 SAMPLES 

load t1b2c5s2.asc 

load tib2c5s3.asc 

load t1b2c5s4.asc 

load t1b2c5s5.asc 

load tib2c5s6.asc 

load t1b2c5s7.asc 

load t1b2c5s8.asc 

load tib2c5s9.asc 

load t1b2c5s0.asc 

accel1=t1b2c5s1; 

accel2=t1b2c5s2; 

accel3=t1b2c5s3; 

accel4=t1b2c5s4; 

accel5=t1b2c5s5; 

accel6=t1b2c5s6; 

accel7=t1b2c5s7; 

accel8=t1b2c5s8; 

accel9=t1b2c5s9; 

accel10=t1b2c5s0; 
% 

accel1=accel1(:); 

accel2=accel2(:); 

accel3=accel3(:); 

accel4=acceil4(:); 

accelS=accel5(:): 

accel6=accel6(:); 

accel7=accel7(:); 

accel8=accel8(:); 

accel9=accel9(:): 

accel10=acceli0(:); 
% 
% 
%THIS SUBROUTINE INPUTS ACCELERATION (g's) and TIME (SECONDS) 
%FROM 2DOF TUNED DECK SIMULATOR TEST BLOWS WITH UP TO TEN SAMPLES 
%INTO THE FUNCTION "SPECTRA.M" IN ORDER TO OBTAIN THE AVERAGE 
%ABSOLUTE ACCELERATION SHOCK SPECTRA. 
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% 
% 
%Output Variable: maxspec=largest value of absolute acceleration computed at each natural 
frequency 
% 
% 
%Calculate the Shock Spectra for Samples 1 through 10. 
% 
cd... 
maxspec1=spectra(accelt,t); 
maxspec2=spectra(accel2,t); 
maxspec3=spectra(accel3,t); 
maxspec4=spectra(accelé.t); 
maxspec5S=sr ctra(accel5,t); 
maxspec6=spectra(accel6,t); 
maxspec7=spectra(accel7,t); 
maxspec8=spectra(accel8,t); 
maxspec9=spectra(accelS,t); 
maxspec10=spectra(accel10,t); 
% 
Freq=[1:1:150]; % FREQUENCY BAND OF INTEREST 
% 
%Calculate the Average Absolute Acceleration Shock Spectra. 
% 
avmaxspec=(maxspec! +maxspec2+maxspec3+maxspec4+maxspec5+maxspec6+ 
maxspec7+maxspec8+maxspec9+maxspect0)/10; 
% 
% 
% 
cd coxfile 
foad t1b2c1s1.ase 
load tib2c1s2.asc 
load t1b2c1s3.asc 
load t1b2c1s4.asc 
load t1b2c1s5.asc 
toad t1b2c1s6.asc 
load t1b2c1s7.asc 
load tt b2c1s8.asc 
load tib2c1s9.asc 
load tib2c1s0.asc 
base1=t1b2c1s1; 
base2=t1b2c1s2; 
base3=t1b2c1s3; 
base4=t1b2c1s4; 
base5=t1b2c1s5; 
base6=t1b2c1s6; 
base7=11b2c1s7; 
base8=t1b2c1s8; 
base9=t1b2c1s9; 
base10=t1b2c1s0; 
% 
base1=base1(:); 
base2=base2(:); 
base3=base3(:): 
base4=base4(:); 
base5=base5(:); 
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base6=base6(:); 
base7=base7(:); 
base8=base8(:); 
base9=baseQ(;); 
base10=base10(:); 
% 
%THIS SUBROUTINE UTILIZES THE FUNCTION "MODEL.M" TO 
%INPUT UP TO 10 SAMPLES OF ACCELEROMETER CHANNEL #1 DATA 
%INTO THE 2DOF TUNED DECK SIMULATOR MODEL TO OBTAIN THE 
%AVERAGE MODELED ABSOLUTE ACCELERATION SHOCK SPECTRA. 
% 
cd. 
[X1ACC1,modspec1]=model(base1 ,t,3,4,500); 
[X1ACC2,modspec2]=model(base2,t,3,4,500); 
[X1ACC3,modspec3]=model(base3,t,3,4,500); 
[X1ACC4, modspec4]=model(base4 ,t,3,4,300); 
[X1ACC5,modspec5]=model(base5,t,3,4,500); 
[X1ACC6,modspec6]=model(base6,t,3,4,500); 
[X1ACC7,modspec7]=model(base7,t,3,4,500); 
[X1ACC8, modspec8]=model(base8,t,3,4,500); 
[X1ACC9,modspec9]=model(base9,t,3,4,500); 
[X1ACC10,modspec1 0]J=model(base10,t,3,4,500); 
% 
%Calculate the Model Average Absolute Acceleration Shock Spectra. 
% 
avmodspec=(modspec1 +modspec2+modspec3+modspec4+modspec5+modspec6+ 
modspec7+modspec8+modspec$+modspect0)/10; 
% 
subplot(311),axis(off) 
text(0.1,1,'°2DOF Tuned Deck Simulator - Test 1') 
text(0.1,0.8,’Spring Beam Configuration: 3 Upper / 4 Lower’) 
text(0.1,0.6,'Test Weight on Deck Simulator: 500 ibs’) 
text(0.1,0.4,'Total Weight on Anvil Table: 3494 Ibs’) 
text(0.1,0.2,,MWSM Hammer Ht: 3 ft (Anvil Travel = 3 in)') 
text(0.1,0.0,'10 of 10 Hammer Drops Used for Average.’) 
subplot(312), plot(Freq,avmaxspec), grid 
titleCAverage Absolute Acceleration Shock Spectra’) 
xlabel(‘Frequency (Hz)') 
ylabel(‘Ave. Accel. (g)’) 
subplot(313), plot(Freq,avmodspec),grid 
title(Modeled Average Absolute Acceleration Shock Spectra’) 
xlabel('Frequency (Hz)') 
ylabel(‘Ave. Accel. (g)') 
set(1,'paperpos',[1.25 2 6 8}) 
print 
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Appendix C.2 (cont.) 


function [maxspec]=spectra(accel.,t) 


NF=~ 50; % NUMBER OF FREQUENCIES 

DF=1; % FREQUENCY INCREMENT 

SF=1; % START FREQUENCY 

Freq=[SF:DF:NF]; % FREQUENCY VECTOR 

wn=2*pi*Freq; 

F=[0;-1]; % STATE SPACE MATRICES 

GG=[1,0]; 

H=[-1]; 

for i=1:NF % STEP THROUGH NATURAL FREQUENCIES. 


E=[0,1;-(wn(i)*2),0]; 

GG=[wn(i)*2,0]; 

[yspec]=lsim(E,F,GG,H,accel,t); 

xspec=yspec + accel; 

maxspec(i)=max(abs(xspec)); 
end 
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Appendix C.2 (cont.) 


function [X1ACC,modspec]=model(base,t,nk1.nk2,_OAD) 
%function [X1 ACC, modspec]=model(base,t,nk1,nk2,LOAD) 
% 
%THIS PROGRAM INPUTS 1 SECOND OF ACCELEROMETER CHANNEL #1 DATA INTO 
%THE 2DOF TUNED DECK SIMULATOR MODEL ONE SAMPLE AT A TIME TO OBTAIN 
%THE AVERAGE MODELED ACCELERATION SHOCK SPECTRA. 
% 
%INPUT VARIABLES: 


% base=anvil table absolute acceleration (base excitation). 
% t=total sampling time period of 1 second duration during which 2048 data points were 
obtained. 


% 

%USER DEFINED VARIABLES: 

% nki=3; 

% nk2=4; 

% LOAD=500; 

% 

% LOAD=weight (500,1000 or 1500 Ibs) atop 2DOF fixture. 

% nk1=number of spring beams (2 or 3) in upper tier of fixture. 
% nk2=number of spring beams (3 or 4) in lower tier of fixture. 
% 

%OUTPUT VARIABLES: X1ACC=model absolute acceleration of mass M1. 


% MODSPEC=model absolute acceleration shock spectra. 
% 
Wt1=738+LOAD; % Weight of Upper Tier in Ibs in 3/4 Configuration 
Wt2=1810; % Weight of Lower Tier in Ibs in 3/4 Configuration 
% Wt1=662+LOAD; % Weight of Upper Tier in Ibs in 2/4 Configuration 
% Wt2=1756; % Weight of Lower Tier in lbs in 2/4 Configuration 
kbeam=2.947e6; % Stiffness per Beam in Ibs/ft 
k1=nk1*kbeam; % Upper Tier Stiffness in Ibs/ft 
k2=nk2*kbeam; % Lower Tier Stiffness in Ibs/ft 
zeta=.035; % Critical Damping Factor 
anvil1=base(: 92. 2; % Accelerometer Channel #1 Data in ft/sec42 
t=t() : % Time in seconds 


% 
%Calculate the mass and mass ratio. 
m1=Wt1/32.2 ; 
m2=Wt2/32.2 ; 
m=[(m1,0:0,m2}; 
alpha=m1t/m2; 
%Calculate Tier Natural Frequencies. 
w1=sqrt(k1/m1); 
w2=sqrt(k2/m2); 
%Calculate System Natural Frequencies. 
var=sqrt((w1%2 + alpha*w1*2 + w242)42 - (2*w1*w2)*2); 
wn1=(1/sqrt(2))*sqrt(w142 + alpha*w1“2 + w242 - var); 
wn2=(1/sqrt(2))*sqrt(w142 + alpha*w142 + w242 + var); 
%Uncouple the equations. 
u11=1; 
u21=(k4 - wn142 * m1)/k1; 
u12=1; 
u22=(k1 - wn242 * m1)/k1; 
U=(u11,u12;u21,u22]: 
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%Determine the Modal Mass Matrix. 
M=U"*m*U; 
%Determine the Modal Stiffness Matrix. 
k=[k1,-k1;-k1,k1+k2]; 
K=U"*k*U; 
%Determine the Force Coefficients. 
FC=(U"[-m1;-m2]); 
%Solve the Uncoupled Equations of Motion. 
% {nddot1}+[2*zeta*wnt]*{ndot1}+[wn142]*{n1}=FC1 1*anvil1/M1 
% {nddot2}+[2*zeta*wn2}*{ndot1}+[wn2°2]*{n2}=FC21 *anvilt /M2 
%o 
%State Space Matrix for Equation #1 
A1=[0,1;-(wn142),-2*zeta*wn1]; 
B1=[0;FC(1:1)/M(1:1)]; 
C1=(1,0;0,1;-(wn1%2),-2*zeta*wn1]; 
D1=[0;0;FC(1:1)/M(1:1)]; 
{y1]=!sim(A1,B1,C1,D1 anvil); 
%State Space Matrix for Equation #2 
A2=[0,1;-(wn2'2),-2*zeta*wn2]; 
B2=[0;FC(2:2)/M(4:4)]; 
C2=[1,0;0,1;-(wn2“2),-2*zeta*wn2]; 
D2=[0;0;FC(2:2)/M(4:4)]; 
[y2]=lsim(A2,B2,C2,D2,anvil1,t); 
% 
%Couple The Equations. 
Y1ACC#ul 1*y1(:,3) + u12*y2¢:,3); 
X1ACC=(11ACC+anvil1)/32.2; % Absolute Acceleration in g's 
% 
%Calculate the Model Shock Spectra for Samples 1 through 10. 
% 


NF=150; % NUMBER OF FREQUENCIES 
DF=1; % FREQUENCY INCREMENT 
SF=1; % START FREQUENCY 
Freq=[SF:DF:NF]; % FREQUENCY VECTOR 
wn=2*pi*Freq; 

F=(0;-1]; % STATE SPACE MATRICES 
GG=[1,0]; 

H=[-1]; 

for i=1:NF 


E=(0,1,-(wn(i)*2),0]; 

GG=[wn(ci)*2,0]; 

[yspec]=lsim(E,F,GG,H,X1ACC,t); 

xspec=yspec+X1ACC; 

modspec(i)=max(abs(xspec)); 
end 
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Appendix C.3: MATHCAD® Worksheet: 


Natural Frequencies 
A._SPRING STIFFNESS: 


System 


The spring stiffness is based on the following beam model: 





Supply the following Dimensions: 


1) Dimensions: Upper Tier Lower Tier 
Length: L} =50 in Ly =50 im 
Load Application Point: a, =13) in ay =13 in 
Area Moment of Inertia: 1, =113 in’ Ty =113 in* 
Young's Modulus: E =30-10° psi E =30-10° psi 


2) Stiffness: 


From the beam bending equation: y ABe (x24 3-a7- 3-L) 


2 


. _ Wa : 
Let x=a and F=W/2: YAR" opi tne} 





K= WwW z 12-E-] Ib 
YAB an(3 L- 4-a) in 





Solving for the stiffness: 


For: Upper Tier Lower Tier 
IZET, ZEIT, 
(Convert inches to feet) KB, =--———____——__l2 KB + 12 


(per beam) 


ay (3Ly- 484) 


KB , =2.947+10° 
ft 
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a2? (Lp 499) 


KB > =2.987-10S 


UpperTier = |= = Lower Tier 


Input Number of Beams Per Tier: Beams ; =3 Beams 5 =4 
K, =KB, Beams, K» =KB 5-Beams 
Tier Stiffness: K , =8.842+10° * K 9 = 1.179+10" . 


B. SYSTEM MASSES: (includes Effective Spring Mass) 





Upper Tier Weight: (W4) Intermediate Weight: (W) Anvil Table Weight: (W3) 
Wy = 1238 W > = 1810 W3 = 4634 
W w W3 
l fal 4 = 2 M 3 = 

32.2 32.2 32.2 

M , =38.447 Ib M5 =56211 bt M3 =143913) hb 

2 2 2 

sec sec sec 
C. TIER NATURAL FREQUENCIES: Upper Tier Lower Tier 

— 
1 [Ki oern lee air? 
| 2-% 4M 1 27 4M > 
F | = 76.326 Hz F 2 = 72.889 Hz 


D. System Natural Frequencies: 


Defining B as: B= 





First Frequency: fy = 61.979 Hz 





Second Frequency: f 5 =115.572 Hz 
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